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Quantum Chromodynamics on a lattice with Wilson fermions and a chirally twisted mass term for two degen- 
erate quark flavours is considered in the framework of chiral perturbation theory. The pion masses and decay 
constants are calculated in next-to-leading order including terms linear in the lattice spacing a. We treat both 
unquenched and partially quenched QCD. We also discuss the phase structure of twisted mass lattice QCD. 



1. Introduction 

Lattice QCD is fighting on several fronts in or- 
der to obtain results relevant to the physics of 
hadrons in the continuum: one wants to con- 
trol the effects (a) of the finite size of the lattice, 
(b) of the non-zero lattice spacing a, and (c) of 
too large quark masses. For small quark masses 
Monte Carlo calculations suffer from a slowing 
down roughly proportional to with p — 2 — i. 
This has so far prevented simulations with light 
Wilson quarks of realistic masses. 

For the extrapolation of the numerical re- 
sults to small values of up- and down-quark 
masses, chiral perturbation theory can be em- 
ployed. It amounts to an expansion around the 
chiral limit at vanishing quark masses. The low 
energy parameters of chiral perturbation theory, 
the Gasser-Leutwyler coefficients, can in turn be 
determined by numerical simulations of lattice 
QCD, see P for a review. 

1.1. Tvifisted mass lattice QCD 

For simplicity we consider QCD with Nf — 2 
flavours of degenerate light quarks, m„ = rud, in 
the following. Consider a "twisted" quark mass 
matrix (T3 acts in flavour space) 

M(w) = m^e'""^^^^ =m + ifij5T3, 

where m = niq cos(u;), fx = uiq sin(ti;). In the con- 
tinuum the angle u can be removed by a chiral 
rotation q = e so that physics does 



not depend on uj. On the lattice, however, chiral 
symmetry is broken even for massless QCD, and 
there is a dependence on to due to lattice arti- 
facts. Twisted mass lattice QCD has been intro- 
duced in [2]. Recently it has been advocated to 
employ a chirally twisted quark mass matrix for 
Wilson fermions in order to improve the efficiency 
of QCD simulations due to full 0{a) improvement 
occurring for tu — Tr/2 [3]. 

1.2. Chiral perturbation theory 

The SU(iV/)L S\J {Nf)ii chiral symmetry 
of massless QCD is spontaneously broken to 
S\J{Nf)v, and in addition explicitly broken by 
non- vanishing quark masses. For Nf = 2 the cor- 
responding Pseudo-Goldstone bosons are the pi- 
ons TTfc. Chiral perturbation theory describes their 
dynamics by means of a low-energy Lagrangian. 
It is expressed in terms of the matrix-valued 
field U{x) = exp(i 7rf,(a;)Tf,/_Fo)j which transforms 
as U ^ LU R^^ under chiral transformations 
{L,R) e S\]{Nf)L®S\J{Nf)R. 

The leading order effective Lagrangian is 

£2 = ^ Tr (a^[/t d^^U) - ^ Tr {xU^ + Ux^) , 

where the quark masses are contained in the ma- 
trix X = 2i?owl. In higher orders, terms with 
more fields and/or derivatives appear, which are 
multiplied by the Gasser-Leutwyler coefficients. 
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1.3. Chiral perturbation theory for lattice 
QCD 

Chiral perturbation theory can also be applied 
to QCD on a lattice, see Oliver Bar's talk at this 
conference. The fields of the low-energy effective 
Lagrangian are still defined in the continuum, but 
the lattice artifacts are taken into account by ad- 
ditional terms in £ proportional to powers of a 
|4l5j . Physical quantities, like and F^, ap- 
pear in double expansions in quark masses (mod- 
ified by logarithms) and in a. In leading order 
the lattice terms read C2L = Tr {p W + 

Up^), where p = 2IVoal, with a new parameter 
Wq. Next-to-leading order calculations have been 
done in 5 up to 0{a) and in 6 up to 0{a?). 

2. Chiral perturbation theory for twisted 
mass lattice QCD 

In view of Monte Carlo calculations in twisted 
mass lattice QCD, it is desirable to extend chiral 
perturbation theory to this case. This has been 
done in (7). The twisting of the mass term x —> 
x{^) — "^Bomq e -"^"^3 can be shifted to the lattice 
term by a chiral rotation p — > p{oj) = e^^'^^p. 

In leading and next-to- leading order (NLO) the 
Lagrangian correspond to the one of [21 with a 
twisted lattice term p{uj). Spurion analysis shows 
that the twisting of the mass term does not pro- 
duce further terms. 

The vacuum corresponds to the minimum of C. 
In contrast to the untwisted case, the minimum 
of the effective action is not located at vanishing 
fields, but at a point where tts 7^ 0, displaying an 
explicit fiavour and parity breaking. Chiral per- 
turbation theory amounts to an expansion around 
this shifted vacuum. 

2.1. Pion mass and decay constant 

Expansion of C in terms of pion fields yields 
the tree level contribution to the pion propagator. 
Loop contributions come from the leading order 
vertices. New vertices arising from the shift of 
the vacuum yield contributions at order a? only. 
As a result we obtain 

™j = Xo + Po + KxVFi){ALl + 2LI - 2LI - LI) 
+8 (xoPo/F2) (4M^- ^ -2Wl~ - 2LI - LI) 



+((xo + po)V(32^2F2)) ln((xo + Po)/A') , 
where xo — 2Bomq , po — 2Woa cos(a;). Here 
and are renormalized chiral parameters and 
A is the renormalization scale. 

At order a, the dependence on the twist angle 
uj shows up in factors cos(a;). It should be noted, 
however, that this is different in higher orders, 
where the above mentioned shift in the pion fields 
introduces new vertices. 

In the case of maximal twist, lu — n/2, the 
lattice artifacts vanish to order a. This has been 
observed for lattice QCD in general in [^, and 
is the basis of the improvement proposal made 
there. 

The second physical quantity we calculated 
is the pion decay constant Ft^, given by 
{0\J'^'°'\TTb{p)) = iF-j^p^^Sab, where J a is the ax- 
ial current. A one-loop calculation gives 

= Fo(l + (4/i;^o')[Xo (2i^ + LI) + po (2Wl 
+W^5'-)]-(l/16^'^o')[(xo + Po)ln((xo + Po)/A2)]). 

3. Partially quenched lattice QCD 

Partially quenched QCD is an algorithmic ap- 
proach to the regime of small quark masses. The 
Monte Carlo updates are being made with sea- 
quarks, which have large enough masses mg in 
order to allow a tolerable simulation speed. On 
the other hand, quark propagators and related 
observables are evaluated with smaller valence- 
quark masses my. Chiral perturbation the- 
ory has been adopted to the case of partially 
quenched QCD in |HE|. 

Combining the approaches mentioned above, 
it appears attractive to simulate QCD with chi- 
rally twisted quark masses in a partially quenched 
manner. For the theoretical analysis of the data 
the extension of the results of [7] to the partially 
quenched case has been done in [TUj . 

In partially quenched chiral perturbation the- 
ory the field U {x) is extended to a graded matrix 
in SU(4|2). The quark mass matrix is in our case 
given by M = diag(mv, my, ms, ms, my, my). 
Its twisted counterpart is M{ijJVtUJs) = 

We have calculated the pion masses mss, "mw, 
mvs and decay constants Fss, Fyv and Fys 
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at NLO (one- loop) including 0(a) TU". At or- 
der a the dependence on the twist angle amounts 
to factors cosujy or coBujg. The expressions can 
be used in the analysis of numerical results from 
Monte Carlo calculations and will aid the extrap- 
olation to small quark masses. 

An extension of the calculation to order a? 
is in progress. For the unquenched case results 
have been published recently by Scorzato and 
Sharpe and Wu [1313] . 

4. Phase structure of twisted mass QCD 

A prerequisite to any numerical simulation 
project is the knowledge of the phase structure of 
the model under consideration. Where are lines 
or points of phase transitions and how do physical 
quantities like particle masses behave near them? 
The phase structure of twisted mass lattice QCD 
has been discussed in 14 11 12 13 on the basis of 
chiral perturbation theory. 

For lattice QCD without twist, Aoki has pro- 
posed the possibility of a phase with sponta- 
neous flavour and parity breaking .15;. An anal- 
ysis of this scenario based on chiral perturba- 
tion theory has been made by Sharpe and Sin- 
gleton A central role plays the potential 
V = — ci Mo + C2 contained in the effective La- 
grangian, where Ug = ^Tr(C/). It contains param- 
eters Ci ^ niq, C2 = 0{a^) for small iriq = 0{a^). 

In twisted mass QCD, where fi = ruq sin(w) ^ 
0, the potential V gets additional contribu- 
tions: V = —ciUq -I- C2 Uq -f C3it3, where C3 — 

Depending on the sign of C2, the possible sce- 
narios are 

C2 > : Aoki scenario near the "critical hop- 
ping parameter" Kc with explicit flavour and par- 
ity breaking, and massive pions. 

C2 < : normal scenario with a 1^^ order phase 
transition extending into the fJ- ^ region, with a 
second order end point at fic = \c2\/F^Bo ^ c? ■ 

On the phase transition line the jump in the 
quark condensate and the neutral pion mass de- 
crease to zero, when the endpoint is approached: 
(A < XX >? - ml, - (l/2F2|c2|)(4ci - ci) ~ 
a\ 

Recent Monte Carlo calculations |16I17| sX (3 — 



5.2 indicate the presence of the normal scenario 
with a first order line. Owing to the associated 
two-phase coexistence and metastability, this rep- 
resents a problem for numerical simulations and 
it would be desirable to have the phase transition 
line as short as possible. 
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Figure 1. Phase structure for the normal scenario 
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